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System and Method For Estimating 
Clock Acceleration and Location Determination 

Jesse Stone 

BACKGROUND OF THE INVENTION 

1. Field of the invention 

This invention relates to a system and a method for processing positioning signals to 
estimate an acceleration of an embedded timing signal, so as to improve the performance of a 
positioning function. 

2. Discussion of the Related Art 

Changes in frequency can reduce the signal-to-noise ratio in a measured positioning 
signal, typically acquired by calculating a correlation between the positioning signal and a 
locally generated replica. This reduction in signal-to-noise ratio is particularly relevant in a 
system in which a coherent correlation is calculated over a long time period. Hence, there is 
a need for a method and a system to estimate the variations of the clock's frequency. In 
addition, there is a need for a method and a system that uses these estimated variations to 
improve processing of the positioning signals. 

The following co-pending U.S. patent applications ("Co-pending Patent 
Applications") assigned to the Assignee of the present patent application disclose relevant 
signal processing techniques: 

1 . Signal Acquisition using Data Bit Information, serial Number 09/888,228, 
filed June 22, 2001, (hereafter, "Application 228"); 

2. Synthesizing Coherent Correlation Sums at One or Multiple Carrier 
Frequencies Using Correlation Sums Calculated at a Coarse Set of Frequencies, 
serial number 09/888,227, filed June 22,2001 (hereafter, "Application 227"); 

3. Extracting Fine-Tuned Estimates from Correlation Functions Evaluated at 
Limited Number of Values, serial number 09/888,338, filed June 22, 2001, (hereafter, 
"Application 338"); 
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4. Determining the Spatio-Temporal and Kinematic Parameters of a Signal 
Receiver and its Clock by Information Fusion, serial number 09/888,229, filed June 
22, 2001 (hereafter, "Application 229"); 

5. Determining Location Information Using Sampled Data Containing Location- 
Determining Signals And Noise, serial number 09/888,337, filed June 22, 2001 
(hereafter, "Application 337"); 

6. System and method to estimate the location of a receiver in a multi-path 
environment, serial number 10/237,556, filed on September 6 2002 (hereafter, 
"Application 556"); and 

7. System and method estimate the location of a receiver, serial number 10/237, 
557, filed on September 6, 2002 (hereafter, "Application 557"). 

The Copending Patent Applications are hereby incorporated by reference in their 
entireties. 

SUMMARY OF THE INVENTION 

The present invention provides a method and a system for a location determination 
based on estimating a clock signal acceleration. According to one embodiment of the 
present invention, the method includes (a) acquiring a first positioning signal; (b) analyzing 
the first positioning signal to provide an estimate of a clock signal acceleration; (c) acquiring 
additional positioning signals based on the estimate of the clock signal acceleration; and (d) 
performing the location determination using the first positioning signal and the additional 
positioning signals. In one embodiment of the present invention, the additional positioning 
signals are acquired using a stacking technique. The first positioning signal may be acquired 
based on a signal-to-noise ratio exceeding a predetermined threshold. 

According to one embodiment, the clock signal acceleration may be estimated by (a) 
dividing the first positioning signal into a predetermined number of segments; (b) estimating 
a phase value for a time point in each of the segments; (c) fitting the phase values into a 
parametric model that depends on the phase values and the clock signal acceleration; and (d) 
deriving the clock signal from the parametric model. Such a parametric model, for example, 
may be based on a constant clock signal acceleration. The parametric model may be based 
on a parabolic function having variables: the clock signal acceleration, an initial phase value 
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and a clock Doppler. In that embodiment, each phase value of the segments may be 
estimated using a quadrature correlation function based on phase values previously estimated. 

According to a second embodiment of the present invention, the clock signal 
acceleration may be estimated by (a) selecting a set of test clock acceleration values and a set 
5 of test clock doppler values based on a clock Doppler estimated from the acquired first 

positioning signal; and (b) choosing one of the test clock acceleration values as the estimate 
of the clock signal acceleration based on evaluating an ambiguity function using the test 
clock acceleration values and the test clock Doppler values. The ambiguity function may be 
based on the magnitude of a complex quadrature phase correlation function. In one 
1 0 embodiment, the test clock acceleration values and the test clock Doppler values are selected 
based on a coarse grid and a finer grid of the clock acceleration values and the clock Doppler 
values. 

The present invention is better understood upon consideration of the detailed 
description below and accompanying drawings. 

15 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a flow chart illustrating a method for acquiring a positioning signal, such 
as a GPS signal, estimating a clock acceleration and, using the acquired signal and the clock 
acceleration, to determining a location, in accordance with one embodiment of the present 
invention. 

20 Figure 2 shows device 200 which is an exemplary circuit suitable for use to carry out 

acquisitions of GPS positioning signals and the location determination described above in 
conjunction with Figure 1 . 

Figure 3 show flow chart 300, which illustrates a first method for estimating a clock 
acceleration, according to one embodiment of the present invention. 

25 Figure 4 show flow chart 400, which illustrates a second method for estimating a 

clock acceleration, according to a second embodiment of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figure 1 is a flow chart illustrating a method for acquiring a positioning signal, such 
as a GPS signal, estimating a clock acceleration and, using the acquired signal and the clock 
30 acceleration, to determining a location, in accordance with one embodiment of the present 

invention. Clock acceleration is a measure of a clock frequency variation occurring while the 
positioning signal, is being captured. The Clock frequency variation is often referred to as a 
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"clock doppler." The Copending Patent Applications (e.g., Application 557) provide some 
examples of ways to perform some of the steps of Figure 1. 

As shown in Figure 1, at step 1 10, a first positioning signal is acquired in any manner 
without using a clock acceleration estimate. The signal-to-noise (SNR) ratio of this first 
5 positioning signal is preferably sufficiently high for reliable acquisition, despite the unknown 
clock acceleration. At step 120, a clock acceleration value is estimated, in accordance with 
one embodiment of the present invention. Such a clock acceleration estimate is illustrated in 
this detailed description by two techniques described in further detail below. At step 130, a 
second positioning signal is selected as a target for acquisition. At step 140, this second 
1 0 positioning signal to be acquired is mixed with a carrier phase term containing a known 
Doppler shift associated with the selected positioning signal, and a varying Doppler shift 
based on the estimated clock acceleration obtained previously at step 120. 

At step 150, if desired, the mixed signal of step 140 may be stacked. At step 160, the 
mixed signal is correlated with a reference code associated with the selected positioning 

15 signal using, for example, a correlation function disclosed in a method described in the 

Application 557. At step 170, the correlation function is compared with a threshold value. A 
signal acquisition is obtained when the correlation function exceeds the threshold value, and 
after one or more validation tests are satisfied. Steps 130-170 may be repeated for several 
positioning signals, until sufficient information for location determination is obtained, or until 

20 all available signal sources have been evaluated. Finally, at step 180, information obtained 
from all the validated acquisitions (e.g., the code-phases of the acquired positioning signals) 
may be used in any suitable manner to obtain a location determination. 

Figure 2 shows device 200, which is an exemplary circuit illustrating acquisitions of 
GPS positioning signals and the location determination described above in conjunction with 

25 Figure 1. As shown in Figure 2, device 200 includes GPS receiver front-end circuit 210, 
which receives a GPS signal via an antenna (not shown). In GPS receiver front-end circuit 
210, the received GPS signal is amplified, filtered, and down-converted to an intermediate 
frequency. Next, the intermediate signal is then sampled and quantized down to two bits 
(e.g., a sign bit and a magnitude bit) at each sampling time point. Conventional GPS receiver 

30 front-end circuits, such as GPS receiver front-end circuit 210, are known to those skilled in 
the art. Each bit of the 2-bit signal samples may be serially shifted into a corresponding one 
of registers 215 to be accumulated as two 16-bit words to allow processor 220. Processor 
220 may write these 16-bit words into a word-addressable memory SDRAM 260 1 via 

1 SDRAM stands for "synchronous dynamic random access memory." In this description, SDRAM 
260 provides a memory system used during program execution. The use of SDRAM integrated 
circuits to implement is merely exemplary. SDRAM 260 can be implemented by any one of many 
memory technologies. 
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memory bus 250. In this embodiment, both GPS receiver front-end 210 and processor 220 
receive a clock signal from oscillator 230. Pulse generator 222 transmits a timing signal to 
processor 220 so as to synchronize device 200 with a host device (not shown). In one 
embodiment, the host device may be a mobile communication terminal, such as a wireless 
5 phone. The host device, in turn, may receive timing information from a network to which it 
is connected. 

In device 200 of Figure 2, host serial interface 224 handles all communication 
between processor 220 and the host device, including, for example, transmitting from the host 
GPS signal acquisition aiding data to processor 220, and transmitting positioning information 

10 back to the host device. Host USB interface 226 and JTAG interface 228 may be provided 
for debugging purposes. Processor 220 processes positioning signals by executing 
instructions of a program stored in Flash RAM 270 2 . Flash RAM 270 also provides non- 
volatile storage for reference signals. The instructions of the program and the reference 
signals in Flash RAM 270 are loaded into SDRAM memory 260 during operation. Memory 

15 SDRAM 260 is also used to store intermediate calculations of the processor 120, such as the 
grid correlation values described in App. 557. 

Figure 3 shows flow chart 300, which illustrates a first method for estimating a clock 
acceleration value, in accordance with one embodiment of the present invention. At step 310, 
a first positioning signal is acquired without using a clock acceleration estimate. The signal- 

20 to-noise (SNR) ratio of this first positioning signal is preferably sufficiently high to allow a 
reliable acquisition, despite the possible presence of an unknown clock acceleration. At step 
320, the acquired signal is divided into n signal segments. For example, in one embodiment, 
a GPS signal of duration 1024 ms is divided into n = 8 segments. Next, at step 330, the IQ 
correlation (i.e., quadrature) values associated with each signal segment of the positioning 

25 signal acquired at step 3 1 0 are calculated. These IQ correlation values each can be used to 
estimate a doppler frequency and a code-phase. If the acquired positioning signal is, for 
example, divided into 8 segments then step 330 generates 8 pairs of IQ correlation values. 
At step 340, phase 0. of each segment is derived from the corresponding IQ correlation 

value: 

30 A _ *tan2(g / ,Z i .) 

W 2k 

This detailed description uses the convention that phase (j>. is associated with the time at the 
midpoint of each segment, denoted (in seconds). As shown in the above equation, the 

2 Flash RAM stands for "flash random access memory." "Flash" refers to a type of non-volatile 
memory device that is electrically programmable and erasable. The use of Flash RAM integrated 
circuits in Flash RAM 270 is merely exemplary. Flash RAM 270 can be implemented by any of a 
number of different non-volatile memory integrated circuits. 
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phase value is expressed in units of cycles, rather than radians, due to the divisor 2n. 
Consequently, with time expressed in seconds, and phase expressed in cycles, the resulting 
estimated clock acceleration is expressed Hz per second. 

At step 350, a parametric model is fitted to obtain the predicted phases 0. . In one 
embodiment, by assuming a constant clock acceleration, a quadratic parametric model is 
used. Fitting to a quadratic parametric model can be carried out in two steps. First, the code 
phases are "unwrapped." This step is performed to resolve the ambiguity in the calculated 
phase of step 340. Unwrapping may be accomplished by the following procedure (for each 
segment, the phase ambiguity is determined based on all previous segments). 

1. For the first segment, the unwrapped phase is accepted as the phase calculated 



2. For the second segment, the unwrapped phase is required to be within one-half 
cycle of the first segment. To obtain the unwrapped phase, an integer number of 
cycles are added to or subtracted from calculated phase; 

3. For the third segment, a straight line fit using the first two segments is used to 
predict the phase of the third segment. The unwrapped phase of the third segment is 
also required to be within 1/2 cycle of the predicted phase; thus, the predicted phase is 
given by: 



4. For the fourth and successively higher segments, the unwrapped phase is obtained 
using unwrapped phases of three previously processed segments. To illustrate, let t 4 
be the time of the segment (i.e., the time at the midpoint of the segment) which phase 
is predicted, t 3 be the time of the segment immediately preceding, t/ be the time of the 
first segment, and t 2 be the time segment halfway between t } and t 3 . Such a selection 
results in a choice of t lt t 2 and t$ spaced apart the maximum amount of time. The 
predicted phase 0 4 satisfies the following (parabolic) equation: 

03 -0 2 (j> 2 -<j) x 04-03 &~02 



above; 



<j> 3 -(j) 2 = (f> 2 -0, 



or 0 3 =(0 2 -0,) 





2 



2 



so that 0 4 is given by: 
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15 



*4 = 



The unwrapped phase is required to be V 2 cycle of the predicted phase. 

After code phases unwrapped, as described above, a parabolic fit is performed on all 
the unwrapped phases. To perform the parabolic fit, the following model phase model is 
5 adopted: 

a,f. <f> 0 are the unknown acceleration, Doppler and carrier phase parameters to be determined 
using the unwrapped code phases. In one embodiment, the parameter values are obtained by 
minimizing the sum of squares of the differences between the unwrapped phases and the 
1 0 phases predicted by the model. In other words, by minimizing: 



20 



Minimization can be accomplished by setting to zero the derivatives of S(a,f,<p o ) with 
respect to a, /and <f> 0 , respectively: 



,2 

To simplify presentation, we define T m = and P m = . Rewriting the above 

equations using T m and P m to obtain: 



\<*T 4 +fr 3 +<j> 0 T 2 =p 2 

\aT l +fT 2 +W x =P x 
±aT 2 +fT i+ <l> 0 n = P 0 



where n is the number of segments. Solving these equations for a,f, <l> 0 , one obtains: 
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(T 4 T 2 -T,T y )(PJ 4 -P 2 T 2 )-(T 4 T x - T 2 T,){P,T 4 - P 2 T,) 
(T 4 T 2 - T,T,)(T 4 n - T 2 T 2 )- (T 4 T X - T 2 T,)(T 4 T X - T 3 T 2 ) 

P x T 4 -P 2 T,-<t> Q (T 4 T x -T,T 2 ) 
(T 4 T X -T 3 T 3 ) 

a = ^(P 2 -fr,-<l> Q T 2 ) 



5 The acceleration value a and Doppler value / can then be used for signal acquisition 

step 140 of Figure 1. 

Figure 4 show flow chart 400, which illustrates a second method for estimating a 
clock acceleration value, according to a second embodiment of the present invention. Steps 
410, 420 and 430 are substantially similar to steps 310, 320 and 330, described above. 

10 Specifically, at step 410, a first positioning signal is acquired without using a clock 

acceleration estimate. The acquired signal preferably has a signal-to-noise (SNR) ratio 
sufficiently high to allow a reliable acquisition, despite the unknown clock acceleration. At 
step 420, the acquired signal is divided into n signal segments. For example, in one 
embodiment, a GPS signal of duration 1024 ms is divided into 8 segments. Next, at step 430, 

1 5 the method obtains in-phase correlation (I) and quadrature phase correlation (Q) values 
associated with each segment of the positioning signal acquired at step 410. 

At step 440, scanning across acceleration space and the Doppler space for a set of 

values for acceleration a and a set of values for doppler /in the vicinity of the measured 

Doppler / of the positioning signal acquired at step 410. Then, at step 450, for each value of 

20 acceleration a (a *) and dopper f(f ') in the selected sets of acceleration and Doppler values, a 
phase value </>. is computed for each segment i: 

a=2w(±a'f, 2 +/7,) 

At set 460, using the computed phase 0,. , an ambiguity function, such as the sum of 

squares of the real and imaginary components of quadrature correlation function, is 
25 computed: 



^(/.cos^+a.sin^) 



-|2 r -»2 



+ 



£(-/ f sin0.+£>.cos^) 
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The pair of acceleration value a ' and Doppler value f that yield the greatest ambiguity 
function is selected to be the predicted acceleration value and Doppler value, respectively, to 
be used in step 140 of Figure 1. 

The number of values of a * and f to be selected at step 440 can be determined by, for 
5 example, the following procedure: 

1 . Define t as the spacing between segments and Tas the total capture 
duration. (For example, if one divides the acquired signal into 8 segments, t would be 
0. 1 28 seconds for a capture duration T of 1 .024 seconds); 

2. For Doppler space (/), determine the number of values using a bin 

10 spacing of 0.25 cycles per captured duration, and search range of 1 cycle per segment 

spacing (e.g., using the t and Tof the previous step, the bin value would be 0.24414 
Hz, and the search range would be 7.8125 Hz, so that 32 values of Doppler f would 
be used); 

3. For the acceleration space (a), a bin spacing of 0.5 cycles per captured 
15 duration squared and search range of 1 cycle per segment spacing squared (e.g., 

using the t and Tof the previous step, the bin value would be 0.47684 Hz/sec, and the 
search range would be 61.035 Hz/sec, so that 128 values of acceleration a ' would be 
used). 

Using the above procedure, for a captured duration of 1 .024 ms, divided into 8 
20 segments, 4096 pairs of (a \f) are used. 

Once the best point acceleration a ' and doppler/ values are identified at step 470 
above, these estimate of a and / may be further refined. It should be recognized that the 
ambiguity function is highly correlated in a and f For example, if the acceleration estimate is 
wrong by an amount Aa, the ambiguity function is maximized by a Doppler / that has an error 

25 of Af= i AaT. One method to address this correlation between a and /is simply to refine 

the search for the peak of the ambiguity function by performing another grid search in the 
vicinity of the best point using a finer grid. For this finer grid, for example 1/8 spacing of the 
previous grid in both the a and / dimensions, a search can be conducted ranges in both the a 
and / dimensions that are each twice the coarse grid values. Thus, the total number of probe 
30 points for this search in the fine grid search is thus 17 x 17 = 289. Alternately, we may use 

derivative-based maximization techniques along with a change of coordinates f=f-^ AaT. 

In this case, the ambiguity function is maximized by searching over values of / and a. 
Because scanning the grid for / and a values may be computationally intensive when the 
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number of segments is large, using eight segments may not be appropriate in instances. In 
general, however, the more segments used, the greater the sensitivity would be, at the 
expense of a longer "time-to-fix" 

This detailed description is provided to illustrate the specific embodiments of the 
present invention described above, and is not intended to be limiting. Many modifications 
and variations are possible within the scope of the present invention. The present invention 
set forth in the accompanying claims. 
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